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Abstract 

Several representations of the exact cdf of the sum of squares of n independent 
identically gamma-distributed random variables X; are given, in particular by a 
series of gamma distribution functions. Using a characterization of the gamma 
distribution by Laha, an expansion of the exact distribution of the sample variance 
is derived by a Taylor series approach with the former distribution as its leading 
term. In particular for integer orders a some further series are provided, including 
a convex combination of gamma distributions for a = 1 and nearly of this type 
for a > 1. Furthermore, some representations of the distribution of the angle <I> 
between (X\,...,X n ) and (1, 1) are given by orthogonal series. All these series 
are based on the same sequence of easily computed moments of cos(<J>). 

AMS 2000 Subject Classifications: 62E15, 62H10 

Keywords: Exact distribution of quadratic forms in non-normal random vari- 
ables, Exact distribution of the sample variance, Gamma distribution, Exponential 
distribution 

1 Introduction 

The distribution of the sample variance s 2 in non-normal cases has attracted sporadic 
attention during the last eight decades. Early investigations, concerning a gamma par- 
ent distribution can be traced at least to Craig (1929) and Pearson (1929), who used 
moment approximations, later investigated more thoroughly by Bowman and Shen- 
ton(1983). For general distributions various approximations were proposed by Box 
(1953), Roy and Tiku (1962), Tan and Wong (1977) and by Mudholkar and Trivedi 
(1981). The latter authors recommend transformations of the Wilson-Hilferty type. 
Exact results seem to be available only for a mixture of two normal distributions (Hyre- 
nius 1950, Mudholkar and Trivedi 1981). See also chapter 4.7 in Mathai and Provost 
(1992). The interest in the distribution of s 2 or 1/s was also stimulated by investiga- 
tions on the estimation of process capability indices, see e.g. Pearn, Kotz and Johnson 
(1992) and Pearn and Kotz (2006). 
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Apart of its intrinsic theoretical value, an analytical representation of the cdf of s 2 
enables more rapid and more accurate evaluations than Monte Carlo methods. Besides, 
it provides a base to investigate the accuracy of the different proposed approximations. 

Now let X\,...,X„ be a random sample with mean X and sample variance 
S 2 = (n — 1 )~ 1 Ya=\ (Xi — X) 2 from a gamma distribution with pdf 

g a (x)=x a - 1 exp(-x)/r(a), (1) 



and cdf 



With 



G a (x)= / gaGM. 



o 



we have 



F:=£X„ Z:=£x, 2 , U 2 :=ZY- 2 , (2) 

i= 1 i= 1 

(«-l)S 2 = {U 2 -l/n)Y 2 , 1/Vn<f/<1. (3) 

By Laha (1954) — see also Lukacs (1955) — the mutual independence of X and the 
sample coefficient of variation S/X = ^/n((U 2 — l/n)/(l — 1/n)) 1 / 2 was shown to 
characterize the family of gamma distributions indexed by a. This is the key to derive 
the distribution of S 2 essentially from the distribution of Z. At first, it follows with the 
angle <t> between (X\,. . . ,X„) and (!,...,!) that 



ff „(r) := Pr{Z<r 2 } = *v(G m (r/U)) 

= &s>(G<xn(rVn cos $)), 
Pi{V^TS<r} = S v {Gan{r(U 2 -l/n)- 1 ' 2 ) 
= #*(G an (rVicot<5)). 



(4) 



(5) 



With the cdf F a „ of U it follows also 

ffa,„(r) = g Y (F a Jr/Y)) 



(6) 



r a " (T(ob) )~ 1 ! n x an - l F a „ 1 
Jo 



Thus 



A« > »W:='- flB, ^(r) (7) 
is the Laplace transform (L.t.) of 

f* ajl (x) := (r(an))- l x a "- l F a „(x- 1 ), < x < y/H. (8) 

A power series for h* a n - and consequently the cdf H a „ - is derived in (fT2] i . . . ( TTol l 
of the follwing section. Due to ©, (O the coefficients of this series determine the 
moments of /* „ and therefore F a/! and the cdf in (0. 
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For a cdf a probability mixture representation is more appealing than a power series. 
Such a mixture of gamma distributions for H a „ is found in d20l in section|2] 

Then the reader might go directly to theorem 14.11 which is the main result con- 
cerning the actual computation of the cdf of S 2 . By means of a few lines of code of a 
computer algebra system some tables of this cdf with at least eight correct digits have 
been computed for some values of a and n. 

Representations for „ by an orthogonal series with Legendre polynomials and 
by a Fourier sinus series are derived too in section [2] Because of the relations (0, (0 
these two series provide two further representations of h* a n orH a n . The corresponding 
series for F a ,„{u) in (|3~T1 i and the orthogonal series derived from (l32l can be used within 
Q. However, the application of theorem l4~Tl is simpler and more accurate. 

In spite of the numerical use of the double series in theorem |4~T1 it would be theo- 
retically more satisfying to have a single alternating power series or even a probability 
mixture representation for the cdf of S 2 . The latter is accomplished exactly for the 
exponential case (a = 1) and nearly for integer a > 1 by theorem |4~2l which contains 
also an alternating power series. However, for numerical purposes, theorem l4~2l is (at 
present) not considered as a competitor to theorem 14.11 since the computation of the 
required coefficients is more cumbersome. The method is explained in section [3] It 
is based on a representation of the cdf of tan<t> by a polynomial on a certain section 
[0, <p„] of the domain of <t>. This formula is found in theorem [3~Tl The proofs of all the 
theorems are given in the appendix. 

Throughout the paper formulas from the handbook of mathematical functions by 
Abramowitz and Stegun are cited by A.S. and their number. The symbol Y,(k) stands 
for summation over all possible decompositions k = k\ + . . . +k„ with non-negative in- 
tegers kj. Moments of positive random variables occur also with non-integer exponents 
and for defective distributions. 

2 The distribution of the sum of squares 

The L.t. of Zi = Xf is given by 

Wait) = (2r(a))-' j°° z a l z -^x V {-^Tz-tz)dz 
Jo 

= (2r(a))-' r((a + ^ )/2) r ("+*)/2 (9) 

= (2r(a))- 1 f -«/ 2 exp(^I) 

( r (?M^-s)->(TW^-s)) 

with Kummer's confluent hypergeometric function M. In particular 

".«>=^»p(iH-£) (io) 
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with the cdf <P(x) = (1 + erf(x/v / 2))/2 of the standard normal distribution. For in- 
teger values (1 — 0[)/2 or (2 — 0t)/2 Kummer's M is given by Hermite polynomials 
(A.S. 13.6.17/18). 

The cdf H a>n (r) = Pr{Z < r 2 } can be obtained by the Fourier inversion formula, 
but some further representations are useful. If y/ a (f) is written as t^ a ^ 2 p a (t 1 ^ 2 ), then 

(vaor=^ a " /2 E/w~* /2 (id 



k=0 



with 



(k)i=l K '- 

= i^Rej^^y- 1 ))"/}^, y = pe i <?. (12) 



p > 0, -n<(p<n, keEo. 

The parameter p was inserted here only for numerical considerations. The j5 a „ £ are 
more quickly computed recursively for « = 1,2,4,8,... than by the integrals. 
Laplace inversion in (fTTT i implies 

By the series expansion of G an in (0]l we obtain with the moments 

»a,n,k := ((VS cos*)""^) - » (a "^ )/2 7a,«i (14) 

ffa ' n(r) = rf^) £ St ^ (15) 

and by comparison with ( TT3b the relation 

( 1) \ fia.n,k 



ifc! ^ ' v ; r((a« + fc)/2)' 



Also the following representations of the distribution functions of £" =1 X, f,S ,U and 
tan<I> are based essentially on the sequence (n a .n,k) or equivalently (y a „ j.)- 

A probability mixture representation for // a .„ by gamma distribution functions is 
obtained as follows: With any scale factor A let be 

lia^,X=^ {m+k) lia,n,k- (17) 

Then, multiplying 

j r \ an - 1 00 

^«'^ A w£ M """" ( - A ' ) " /t! 
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by e- lr ZJ =0 (Xr) e /£i we obtain 

Xg an {Xr)j^8 a ^ x {Xrf/k\ (18) 
with the alternating k th order differences 

8 a ,n,k,X = (-A)Va,»,0,A = £ Q (19) 



and consequently 



f CCn + k- 1\ 

#a,n( r ) = L <W,*,A ^ J G an+HAr). (20) 



In particular with A = ^fn this is a probability mixture since 
S a ,n,k^n = ^((cos4>)«"(l - cos4>)*) > and 



E( w ft * -1> )w=i (2D 



(t=0 



by virtue of the binomial series with < 1 — cos<f> < 1 — 1/ \fn. However, for numerical 
approximations different values of A should be more suitable, e.g. 



(2r(a))« r(an/2) 



with jU a „ a = 1- 



Besides, two further representations of Z/ a „ are obtained from the following or- 
thogonal series for /*„ from 0, ©. With 



f k 

/ xfa„(x)dx = 
JO 



T(an) an + k' 
implied by (TT4l) . ( TT5l >. and the shifted Legendre polynomials 

k r \ 

(A.S. 22.2.11), it follows 



(23) 



P*k(y) = LK,y\ 0<y < 1, J\p* k {y)fdy= , (24 ) 



V « i =0 



(25) 



c a,n.k 



an/2 £ v 

1 >V{an)p Pk Jan + j 
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Thus, with (A.S. 1 1.4.26) and the modified spherical Bessel functions (A.S. 10.2.2) we 
find 



H afl {r) = r a "exp {~\r^ £ (-l) k c a , n , k ^ I k+l/2 {^ r ^) ■ (26) 



1/2 

-l)*e OA * ( Tlz ) 

k=0 

Consequently, the Fourier transform of f* a n is representable by 



= ex P i^^*) E [kc a,n.kk (27) 

with the spherical Bessel functions from (A.S. 10.1). For absolutely large real t this 
series is numerically more suitable than the power series from dl5) . Such values are 
required as coefficients in the following Fourier sinus expansion for /* n . 



— - ) , < x < , (28) 

.... v 

with 



ba, n , n = -^Im ("^) 



(29) 



2 ^ (-l)(-+*-i)/2 Cai „ >ftA(m ^/2) 

m+A: odd 



Because of lemma IAT21 the function has a square integrable derivative at least for 
min(an, n) > 2 which entails the absolute uniform convergence of the orthogonal series 
in d25l > and d28b . Inserting d28l into © leads to 



H a , n {r) = r a " £ ™ UV " 2 b a ,n, m (l - (-l) m e- r ^) . (30) 



= i" 

Finally, with the functions 

sf ' (C):= £r(/i+*/2) 

= (r(^))- 1 M(l,/3;C 2 ) + (T(/3 + l/2))- 1 CM(l,^ + 1/2;C 2 ), P > 0, 
= C 2(1 ^>/3-i(CVG/5-i/2(C 2 )), if P>h 
the following integral representation can be derived from (ITTb . ([T2] >, (fOt : 

*Mr) = — I Re{(j3«(y- 1 )) Sfi +eD ,/ 2 (ry)}«ty, 
y = pe'f, p>0. 



6 



3 The distribution of the angle <3> 



Two representations of the cdf F a<n of U — (-^/ncosO) 1 follow directly from ( fZSb and 



F a , n {u) = r(an)u a "- 1 ^ r Y,c aM P*(l/(V^u)) 

(31) 

mK 



Y{an)u an V fc a ,„, m sin —= 



Besides, an orthogonal expansion with Legendre polynomials for the density of U 2 is 
obtained by the moments 

^ ) -r(an+2fe)^I = lT(^r' (32) 

derived from @. 

Direct use of these orthogonal series within (f5]) is possible to get the cdf of S, 
but the double series in theorem [4. II is numerically more favourable. As mentioned 
in the introduction it would be theoretically more satisfying to have a single alternat- 
ing series or a probability mixture for this cdf. A single power series is not directly 
available by a power series expansion of G an in © since the moments of cot<t> do 
not all exist. However, at least for integer a, this obstacle is avoided by splitting 
the domain [0,arctan(v« — 1)] of O by <p„ = arctan(l js/n— 1), which is the max- 
imal angle <p for which the whole cone {<t> < <p} is contained completely within 
{x\ , . . . , x„ \x\ , . . . ,x„ > 0}. For integer a the theorem lXTI below provides a polynomial 
representation of the cdf of tan <t> restricted to < tan<p < (« — l) -1 / 2 . This simple rep- 
resentation is used within (01 to integrate over tan <p < (n— l) -1 / 2 . The integration over 
tan<p > (n — l) -1 / 2 needs only truncated moments of cot<I> = (1 — cos 2 <I>) _1 / 2 cos<l>, 
obtained by the binomial expansion with the truncated moments of cos <t>, given below 



Theorem 3.1 Let a be a positive integer, then 
W a ,„(r) := Pr{tan4><f} 



(33) 



nan) [(a -£ n ' 2] f"- 1+2 > „ , „ 1/2 

E a a ,„,2j — r~^ ,0<t <(n-l) x l\ 



(r(a))"n««/ 2 ; to ' - 1 n-\+2j 
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where the coefficients aa,n,2j are given as the unique solutions of the linear equations 



[(«-l)«/2] 



j=0 



£ r[-+pn+m-j)r 



n-\ 



J a a,n,2j : 



(34) 



f 2Ml£nr(l/2 + J 3+m y )/(2m y )! 

(m);=l 



a = 2/3 + 1 



2 m+ J' i:ri r ( 1 /2 + ^ + ^)/(2m, + l)!, a = 2/3 

" (m)./=l 

m,; = 0,...,[(a-l)n/2]. 
With any p > the sums within the right hand sides of d34l i are also given by 



1 



p n n Jo 



a 



1 



P'"7t Jo 



Re 



2 

a + l 



Q -nntp 



M 



a + l 3 1 



2 '2 ; 4 Pe 



d<p, a = 2j3 + l, 



d<p, a = 2)3, 



(35) 



After Kummer's transformation (A.S. 13.1.27) the above confluent hypergeometric 
functions M can also be expressed by Hermite polynomials due to (A.S. 13.6.17/18). 

The truncated moments fa.n,k = S'((cos<P) an+k ), defined by integration over 
<P > <Pn, follow from theorem IXTl bv straightforward calculation. 



Y{an) 1 ' ( "-^"/ 2 ' 

Ja.n.k - Ja.n.k (Y( a ^n n an/2 2 f- fl ° ! '"' 2 - /Z 



n-1 , ■ (g-l)n+fr+l _ •. 1 
2 2 ■/'« 



(36) 



with the incomplete beta function B(a,b;x) — f^t" l (l—t) b l dt. In particular with 
a = 1 we have 



{n-iy. n^-^l 2 (n-\ k+l 1 



(37) 



Two further remarks: The equations (f34t can also be solved by explicit matrix in- 
version. The matrix (r(A +j3« + m — j*)), m,j = 0,...,N = [(a — l)«/2], has the 
structure T(jc)Diag( . . . , (x) m , . . .) ((x + m)N-j) with factorials (x) m — x(x + 1) . . . (x + 
m—l). The inverse (p m j(x)) of ((jr+m)#_y) contains only polynomials p m; - of degree 
m in its m th row. These polynomials can be obtained by interpolation from the m+ 1 
values 



(-!)/+*+! 



Pmj {-{N-k)) 
(k = Q,...,m) 



(j + k-m)\(N-(j + k))\ 




m—k<j<N—k 
otherwise 



> . 
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The cdf F\ M of Ui,„ can be applied to a goodness of fit test of the hypothesis Hq : F = Fo 
against H\ . F 7^ F), where Fo is any specified continuous cdf. Under Hq the test statistic 
T!j=i (Fo(Xi:k) — FoQ^i-i:k)) 2 > obtained from an ordered sample X^, i— l,...,k = n—l, 
(Xo-.k = X n ± = oo), has the same distribution as Uf n . The asymptotic normality of 
U 2 n can be derived from the asymptotic bivariate normal distribution of ((T„ — «)/ y/n, 

(Z„ — In) I \/2n) with correlation p = 2/y/5. The convergence to a normal distribution 
is slow and might be accelerated by a suitable transformation. 



(38) 



4 The distribution of the sample variance 

With the functions 

H a ,n,m(r) ■= £{Gan+ m {r\/n cos®)) 

an+m oo ,. I \k 

r(an + m) j^an + m + k k\ ' 

m e N , H a ,„. k = <f (V«cos <j)) a " +k = n (a " +A)/2 7a,„,A- from (O we obtain 

Theorem 4.1 Let X\,...,X n be independent random variables with den- 
sity (T((x))~ l e~ x x a ~ l , a > 0, then the cdf of the sample variance S 2 — 
{n-\)- l i:U(X-X) 2 isgivenby 

=a^_ Jl _f T(a« + 2;) ( d V' 



Pr{(«-l)^<-^} = E^P^j 09) 
£ f £ ^(an + fc)/2 + J - l^j (-r)* ^ 1 (4Q) 



r(an)jto\^ V ./' J an + k k\ 

t (t W,A ( an \ k ~ 1 ) f -V (41) 



y=0 \i:=0 V 

vv/f/i f/ze alternating k th order differences 

X - V r i y + ^)/2 + J - A Ma,n,2j+l , 49 . 

l W/i.A - Lt _1 M £ ) I y ) A a„+2j+/ (42) 



/ = 1 



and any X > 0. 77ie series in ( 1591 ) 25 absolutely convergent with terms of the order 
O (j~( n+ ^/ 2 ^j, where the O-constant depends only on a andn. 

As a simple numerical example we obtain with a = 1 and n = 10 the value 
Pr{S < 2} = 0.98530379 . . . with at least eight correct digits. 

Finally, the method, explained before theorem [3~Tl is applied to obtain the alterna- 
tive representations in theorem l4~2l for integer a, including in particular a probability 
mixture for the exponential case a = 1 . The degree of the polynomials in theorem [3~TI 
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and the rate of convergence of the binomial series in (l43t limit the numerical use of 
theorem l4~2l to moderate values of an. 
The truncated moments 

M a , n , k = i ((«*•)«**) = t ( (a " + k) f + ? a , n , +2] (43) 



are used with y a . n ,k from d36l >. obtained by integration over q>> (p n = arccot(\/« — T). 



Theorem 4.2 Wi'f/i f/ie truncated moments from (|45J) ana! ( 1561 ), f/ie pdfw a ,n and the cdf 
W a ,n o/tan<I> on < <p < q> n frotn theorem \3.1\ the distribution of S is given by each of 
the following three formulas (O, (O, d?6l ): 

Pr|\/«-15<r| 

_ (44) 

with the polynomials 

/wM=E V(a " +t,/2 ^ 2 ' 



;=0 



(,-»"" y M a , n . k (- ry /n) k 
F(an) ^ Q an + k k\ 

+ W a Jl/Vn^T)G a „ (ry/n(n-l)) + / g a „ (jc)Wa- (r^/n/x) dx. 

V / JrJn(n-\) 



The infinite series in ( |45| ) is a/so given fey 



LV^^Y ^ G an+ jt(Ar) (46) 



wifn any X > ana" fne alternating k tb order differences 



K.n.k.X = ("A)*Ma^ = £(-!)' (J) 



1=0 



Remarks: In particular, with a = 1 the only coefficient in wi „(f) is ai „o = 
(n - l)b n -i with the volume 6„_i = K^^/T (2±i) of the (n - l)-unit ball. Then 
the last two terms in d45T > are reduced to 

Mn£-W^)^ G ^'S^l (4?) 
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The truncated moments fln,k were given in d37b . Only positive coefficients & an kX 
arise in d46]> with A = \Jn{n - 1) since 1 > 1 - -3= > 0. The last term in (@3]l is 



bounded by W a ,„(l/\/n— 1) (l — G an {ry/n{n — 1))J, which is often neglectible. For 
actual computations the choice 



A = Jn 



<#((cot<I>) a ") ^ 1/(an) 



l-Wa,„(l/Vn^T) 



might be more favourable. With this X the leading term in ( f46b becomes 
(1-W a ,„(l/V^ T))G an (Ar). 

It would be desirable to find a similar probability mixture representation for the cdf 
of S 2 for all a > but presumedly, also in case of its existence, the computation of the 
required coefficients would not be easy. 

Acknowledgement: The author would like to thank S. Kotz for his strong interest 
in this work and some hints at the literature within a personal communication. 



Appendix 

Proof of theorem UJ] Let a be any positive integer and CO the Lebesgue measure on 
the (n — 1 )-unit sphere S„_ i . Then 

ce-1 



( \ 



(A.l) 



2T{{an+k)/2yh\YX\^{a + kj)/2)/kj\ , a + it; odd. 
Wi =1 



After an orthogonal transformations = Ty with y\ = rT x l 1 YIi=\ x i — cost?i, followed 
by transformation to polar coordinates, we obtain with c; = cos r?y, Sj = sin the same 
integral value as 

n 



= n k l 2 



s„-i i =0 ;=i V *-./' / / 
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By comparison with dA.U it follows 



[(«-J>/ 2 l /( _i) B + j fc + i \ fn-l , 



.7=0 



(A.3) 



:2«-*/ 2 &!£n r ((«+W2)A;' 



With a = 2j3 + 1 only even kj = 2m j occur with sum k = 2m, and with a = 2j3 only 
/fy = 2m / + 1 with sum k = 2m + n, which leads to the linear equations in (|34| |. 
If I is odd then the coefficients a a „ c vanish since 

is a linear combination of the integrals 

»-2 \)=2 / .7=2 

which are different from zero only for even m2,...,m„, Y!j=2 m j = L 1= i^j = ^- 
If f = tan<p < (n - l) -1 we obtain with S„_i.<p = S„_i n {#i < <p} that 



Pr{tan4><f} = (r(a)) _n / We^'x^dxj 

M<<P}y=i 



(T(a))- n / / p^-'expC-pv^cOn^i/V^ + ^iT,)"- 1 ^ </p 

1 l^U«; V „ / n-2+2; -n-2y 



(r(a))"« a "/2 



a,n,2; / s \ Cj W T>i , 



where the integrals are given by (tan<p)" 1+2/ /(n — 1 +2 /', which provides the asserted 
result in □ 

Comparing the series 

~ (_!)* (z/y) («n^)/2 



with the corresponding one for T(an + 2 j)G an+ 2j(\/z/y) we obtain the identity 

-#-V G an = f #Y G an+2 , fv^) (A.4) 



<9y/ \ J Y{an) \dz 
by differentiation. 
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Lemma A.l 



dy 



fzfy 



<— — — , y,z>o. 

% jyJ 



(A.5) 



Proof. Using (IA.4t and Cauchy's integral formula for the derivatives of 

^T(an + 2j)G an+2j (V~z) = l -z an ' 2+] - 1 exp(-Vi) 
the bound in JA.5t follows from 



(;-l)! 



r(an+2j)G an+2 ;(Vz) 



= lim 



< — J y a "/ 2 - 1 exp(- v ^72) £ /y = 2 an / 2 r(an)/7r 
where the way of integration is 

w R = {Re ,,p | - | < <p < |} U {iy|* > y > -fl} . 



□ 



Lemma A.2 Lef 4> denote the angle between (X\, . . . ,X„) and (1, ... ,1) from (0, 
W a ,n(t) the c4fo/tan4> and b n -\ = ^'" _1 '/ 2 /r(^) f/ze volume of the (n — \)-unit 
ball, then 



W a Jt) 



T{an)b n ^ 



(r(a))"n( a - 1 /2)» 



one/ m particular 



Furthermore 



W hn (t) = ^—^i bn -,t n -\ 0<t< (n- l)-'/ 2 . 



(A.6) 



(A.7) 



(A.8) 



Proof. Let x,x g denote the arithmetic and the geometric mean of {x\,... ,x„ > and 
<p the angle between (xi, . . . ,x n ) and (1, . . . , 1). Setting Xi = (1 + e,)x with any fixed x 
we have the relation 



x 2 ( x i ~ x ) 2 = J! e ? = n tan2 *P 



and for <p J. 



- x-qo + e <) ~ x " ex p ( 5 X>n ~ 
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It follows 

W a ,n{t) = Pr{* <(?}=( UgaiXiidXi 

J<S><<pf = l 

= (r(a)) - " / x"/ a_1 'exp(— nx)dx\ . . .dx n 
Jt&<(p 

Jo 

which provides ( IA.6I ). 

For a = 1 this asymptotic relation can be replaced by an equation if <p < q>„ = 
arctan(l /y/n — l), which is the largest angle <p for which the whole cone {<t> < <p} is 
contained within {x\ , . . . ,x„ \x\ , . . . ,x„ > 0}. 

To prove (I A.81 > we obtain from ( IA.6t 

Pr{|ln(cos4>)| < e} = Pr{ln(l + tan 2 4>) < 2e} 

~ W an (V2i)~c a ,„e (n - i y 2 , e|0, 

with a factor c a ,n determined by ( IA.6t . For 5 — > °° the relation (IA.8b follows from 
the Hardy-Littlewood-Karamata Tauber theorem (cf. e.g. chapter 13 in Feller 1971), 
applied to 

<f((cos$)' s ) =S (eq»(-sr|ln(cos*)|)). □ 

Proof of theoremWJ\ With any G (0,1) #t/ (G a „((z/(£/ 2 - ?V»)) 1/2 )) can be 
represented by the Taylor series 

s>((-i)^(^)u)(?)' 

With I/ -1 = -y/ncosO and the bounds from Lemma |A. 1 1 the terms of this series are 
absolutely bounded by 

-yan/2 

g ((cos4>) 2i ) &/j = 0(^r ( " +1)/2 ) 

K 

due to d A.81 > with s — 2j and on (9-constant depending only on a and n. 
For $ t 1 it follows from © and ( IA.4I ) that 

Pr{(„-l) S *<z = ,*} = ft (g„ (,/=?=)) 

- E^*((|)'Wv^)) 

p T{an)j\nJ \dzj 
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with the functions H a „2j from 
With any A > we obtain 



d ( T{an + 2j) ( d V 



H a ,n,2j(Vz) 



dr \ T(an)j\ni \dz, 

^ ( (an + k)/2 + j- l\ jlg.n.lj+k (-1) n ^an-l+fe/ 5 2/ y 

LI / Ja«»+2;+* it! (Ar) ^ • 



After multiplication by exp(— Ar)^^° =0 (Ar)^/£! the series d4lT > is obtained by integra- 
tion over r. □ 



Proof of theorem WM With the pdf w a , n and the cdf W a ,„ of tan<t> from theorem l3TTI 
and the truncated moments M a .,i,k = & ((cot<fr) an+k ) from d43l it follows with © that 

PrjV^ - TS<r} = g (G„„(rVn cot*)) 



/ Gan{ryfc/t)wa*{f)dt+ \ £ 

Jo r an S, 



/o ' r(ara) j~ Q an + k k\ 

After integration by parts and the substitution f — > r^fnjx the above integral becomes 

W a „ (l/Vn - 1) G a „ (ry/n(n-l)) + / gan(x)W an (ry/n/x) dx. 

Inserting the series for M a n ^ from d43l formula (l44l > follows by rearranging the 
resulting absolutely convergent double series according to the truncated moments 
fa,n,k = § ((cos<t>) an+k ) from @5J. 

The series in (|46*T > is obtained by the same way as (|2T1 i. □ 
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